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Abstract 

We tabulate values of the radial Schrodinger wave function or its first nonva- 

nishing derivative at zero quark-antiquark separation, for cc, cb, and bb levels 

that lie below, or just above, flavor threshold. These quantities are essential 

inputs for evaluating production cross sections for quarkonium states. 
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Fragmentation of partons into quarkonium states recently has been recognized as an 
important component of quarkonium production in high-energy collisions Thorough 
investigation of this new source, and others, is made timely by the CDF Collaboration's 
observation that conventional sources substantially underestimate the yield of prompt J/if> 
and if)' in 1.8-TeV pp collisions [0. In the case of the if>', which is not fed by the cascade 
decay of the narrow \c states, the observed cross section exceeds the calculated one by more 
than an order of magnitude. 

Calculation of the production rate by fragmentation can be separated into a parton-level 
piece that can be evaluated using perturbative techniques and a hadronic piece expressed in 
terms of the quarkonium wave function. We have earlier tabulated the values at the origin of 
the radial wave function, or its first nonvanishing derivative, for narrow levels of the be (B c ) 
system ||. These have been used to estimate the B c production rate j4j. Here we present 
the corresponding information for the cc {J/if>) and bb (T) families. Although many of these 
numbers have appeared in the literature when quarkonium spectroscopy was in flower, they 
usually were given implicitly in calculations of leptonic widths or similar observables. Our 
purpose here is to record the relevant properties of the wave functions of all the narrow 
levels in a form convenient for evaluating cross sections for quarkonium production. 

We consider four functional forms for the potential that give reasonable accounts of the 
cc and bb spectra: the QCD-motivated potential given by Buchmiiller and Tye f5[, with 

m c = 1.48 GeV/c 2 m b = 4.88 GeV/c 2 ; (1) 

a power-law potential ||, 

V(r) = -8.064 GeV + (6.898 GeV)(r • 1 GeV) ai , (2) 

with 

m c = 1.8GeV/c 2 m b = 5.174 GeV/c 2 ; (3) 
a logarithmic potential |7j, 
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V(r) = -0.6635 GeV + (0.733 GeV) log (r • 1 GeV) , (4) 

with 

m c = 1.5 GeV/c 2 m b = 4.906 GeV/c 2 ; (5) 
and a Coulomb-plus-linear potential (the "Cornell potential") M, 

V(r) = ~ + L , (6) 



r a 2 



with 

m c = 1.84 GeV/c 2 m b = 5.18 GeV/c 2 (7) 
K = 0.52 a = 2.34 GeV" 1 . (8) 

For quarks bound in a central potential, it is convenient to separate the Schrodinger 
wave function into radial and angular pieces as \l/ n ^ m (r) = R n i(r)Y em (9, 0), where n is the 
principal quantum number, I and m are the orbital angular momentum and its projection, 
R n e{r) is the radial wave function, and Ye m (6, 4>) is a spherical harmonic 0. The Schrodinger 
wave function is normalized, / <i 3 r|\I/ n £ m (r)| 2 = 1, so that J °° r 2 dr\R ni {r)\ = 1. The value of 
the radial wave function, or its first nonvanishing derivative at the origin, 



(9) 



r=0 

is required to evaluate production rates through parton fragmentation. The quantity 
|i?^(0)| 2 is presented for four potentials in Table | for the narrow charmonium levels and 
in Table |I| for the narrow Upsilon levels. For ease of reference, we reproduce in Table [1 1 1 



our predictions for the B c wave functions, with some computational improvements in the 
entries for the Cornell potential. The strong Coulomb singularity of the Cornell potential is 
reflected in spatially smaller states. 



In view of the efforts ]I0| to resolve the ip' anomaly with cascades from above-threshold 
states, we quote values for the cc 3D, 3P, and 3S levels that lie near 3.8, 3.9, and 4.0 GeV/c 2 , 
respectively, and for the bb 4F, 4D, 4P, and 4S levels that lie near 10.35, 10.45, 10.52, and 



10.6 GeV/c 2 , respectively. It is likely that these will be significantly modified by coupled- 
channel effects gin]]. 

If all the potentials describe the cc and bb observables, what accounts for the wide 
variation in the values of |i? n o(0)| 2 and the corresponding quantities for orbital excitations? 
The leptonic decay rate of a neutral (QQ) vector meson V° is related to the Schrodinger 



wave function through [PL2| , pL3 1 



where N c = 3 is the number of quark colors, eg is the heavy-quark charge, and My is the mass 
of the vector meson. The QCD correction reduces the magnitudes significantly; the amount 
of this reduction is somewhat uncertain, because the first term in the perturbation expansion 



is large ||14|| . Each of the potentials we use corresponds to a particular interpretation of the 
radiative correction to the leptonic width. Similar effects may enter the connection between 
wave functions and fragmentation probabilities. 
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TABLES 



TABLE I. Radial wave functions at the origin and related quantities for cc mesons. 



Level 




\R { >)\ 2 








QCD (B-T), Ref. |] 


Power-law, Ref. || 


Logarithmic, Ref. 


Cornell, Ref. § 


IS 


0.810 GeV 3 


0.999 GeV 3 


0.815 GeV 3 


1.454 GeV 3 


2P 


0.075 GeV 5 


0.125 GeV 5 


0.078 GeV 5 


0.131 GeV 5 


2S 


0.529 GeV 3 


0.559 GeV 3 


0.418 GeV 3 


0.927 GeV 3 


3D 


0.015 GeV 7 


0.026 GeV 7 


0.012 GeV 7 


0.031 GeV 7 


3P 


0.102 GeV 5 


0.131 GeV 5 


0.076 GeV 5 


0.186 GeV 5 


3S 


0.455 GeV 3 


0.410 GeV 3 


0.286 GeV 3 


0.791 GeV 3 
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TABLE II. Radial wave 


functions at the origin 


and related quantities for bb mesons. 


Level 




l<i(o)l 2 








QCD (B-T), Ref. §] 


Power-law, Ref. || 


Logarithmic, Ref. 


Cornell, Ref. § 


IS 


6.477 GeV 3 


4.591 GeV 3 


4.916 GeV 3 


14.05 GeV 3 


2P 


1.417 GeV 5 


1.572 GeV 5 


1.535 GeV 5 


2.067 GeV 5 


2S 


3.234 GeV 3 


2.571 GeV 3 


2.532 GeV 3 


5.668 GeV 3 


3D 


0.637 GeV 7 


0.892 GeV 7 


0.765 GeV 7 


0.860 GeV 7 


3P 


1.653 GeV 5 


1.660 GeV 5 


1.513 GeV 5 


2.440 GeV 5 


3S 


2.474 GeV 3 


1.858 GeV 3 


1.736 GeV 3 


4.271 GeV 3 


4F 


0.414 GeV 9 


0.627 GeV 9 


0.456 GeV 9 


0.563 GeV 9 


4D 


1.191 GeV 7 


1.396 GeV 7 


1.119 GeV 7 


1.636 GeV 7 


4P 


1.794 GeV 5 


1.593 GeV 5 


1.377 GeV 5 


2.700 GeV 5 


4S 


2.146 GeV 3 


1.471 GeV 3 


1.324 GeV 3 


3.663 GeV 3 


5F 


1.075 GeV 9 


1.302 GeV 9 


n hq/I c^awQ 
u.o»4 (jev 


1 ^90 PoV9 
l.OZU 


5D 


1.722 GeV 7 


1.689 GeV 7 


1.289 GeV 7 


2.417 GeV 7 


5P 


1.935 GeV 5 


1.504 GeV 5 


1.252 GeV 5 


2.917 GeV 5 


5S 


1.956 GeV 3 


1.231 GeV 3 


1.077 GeV 3 


3.319 GeV 3 



s 



TABLE III. Radial wave functions at the origin and related quantities for cb mesons. 



Level 




\R>)\ 2 








QCD (B-T), Ref. §] 


Power-law, Ref. 


Logarithmic, Ref. 


Cornell, Ref. § 


IS 


1.642 GeV 3 


1.710 GeV 3 


1.508 GeV 3 


3.184 GeV 3 


2P 


0.201 GeV 5 


0.327 GeV 5 


0.239 GeV 5 


0.342 GeV 5 


2S 


0.983 GeV 3 


0.950 GeV 3 


0.770 GeV 3 


1.764 GeV 3 


3D 


0.055 GeV 7 


0.101 GeV 7 


0.055 GeV 7 


0.102 GeV 7 


3P 


0.264 GeV 5 


0.352 GeV 5 


0.239 GeV 5 


0.461 GeV 5 


3S 


0.817 GeV 3 


0.680 GeV 3 


0.563 GeV 3 


1.444 GeV 3 
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